multipleOlogit <- function(formulas,y=data,modelnames,...){
	require(MASS)
	foo <- vector("list",length(formulas))
	for(i in 1:length(formulas)) {
		foo[[i]] <- polr(substitute(i,list(i=formulas[i])), 
                     data=y,Hess=TRUE)
		foo[[i]] <- summary(foo[[i]])[1]$coefficients
		foo[[i]] <- cbind(foo[[i]],
						foo[[i]][,1]-1.96*foo[[i]][,2],
						foo[[i]][,1]+1.96*foo[[i]][,2])
		foo[[i]] <- data.frame(foo[[i]],modelnames[i],rownames(foo[[i]]))
		colnames(foo[[i]]) <- c("Estimate","SE","T-Score",
                            "Y0","Y1","Model","Var")
			}
	return(foo)
}
